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Abstract 

This paper compares two methods of statistical mechanics used to study a classical 
Coulomb system 5" near an ideal conductor C. The first method consists in neglecting the 
thermal fluctuations in the conductor C and constrains the electric potential to be constant 
on it. In the second method the conductor C is considered as a conducting Coulomb system 
the charge correlation length of which goes to zero. It has been noticed in the past, in 
particular cases, that the two methods yield the same results for the particle densities and 
correlations in S. It is shown that this is true in general for the quantities which depend 
only on the degrees of freedom of S, but that some other quantities, especially the electric 
potential correlations and the stress tensor, are different in the two approaches. In spite 
of this the two methods give the same electric forces exerted on S. 
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1. Introduction 

In the equilibrium statistical mechanics of classical Coulomb systems (for instance 
electrolytes), sometimes one is led to assume that some wall is an ideal conductor (for 
instance for mimicking an electrode) . Two methods have been used to deal with a classical 
(i.e. non quantum) Coulomb system S near an ideal conductor C. The first one is to 
consider from the beginning that the conductor C is ideal and take this into account by 
constraining the electric potential to be constant on C [1-3] . The second method is to treat 
the conductor C as a genuine Coulomb system with a microscopic structure and take the 
limit of zero correlation length [4, 5] ; indeed, in that limit, the statistical average of the 
charge density on the conductor C becomes a surface charge density of zero thickness, a 
characteristic feature of ideal conductors. 

Both methods gave the same results for some quantities in the Coulomb system S, 
such as the particle densities, and even the fiuctuations of these densities as described by 
the particle correlation functions. The reason for this agreement about some fluctuations is 
not obvious. In the first approach, there are no fiuctuations inside the conductor C. In the 
second approach, there are thermal fiuctuations inside the conductor C ; for instance, the 
potential-potential correlation function has a universal simple form [6] (in 3 dimensions, 
ksT divided by the distance) for distances large compared to the microscopic scale and 
these correlations do not disappear as the zero correlation length limit is taken. However, 
this seems to have no infiuence on some quantities in the Coulomb system S. Neverthe- 
less we shall show that some other quantities in S' - for instance the electric potential 
correlations - are different depending on the method used to compute them. 

The aim of the present paper is to discuss the relationship between the two approaches. 
We first treat a simple two-dimensional example in section 2, where we consider two parallel 
lines ; one is the ideal conductor C and the other one is the Coulomb system S. The two 
methods are worked out to find the electric potential and charge correlations and the 
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stress tensor, and we compare the results. In section 3, we treat the general case, look at 
quantities such as the partition function, the correlations and the stress tensor, and again 
we compare the results obtained through the two methods. 

The reader may either look first at section 2, or go directly to section 3. 

2. A simple two-dimensional example 

2.1. The model 



For simplicity, we consider a system of "restricted dimension" [7]. The Coulomb 
interaction between two point-charges q and q' at a distance r from each other has the two- 
dimensional form —qq'ln r, but the particles are constrained to live on one- dimensional 
lines. In the plane xOy, the line y = is the conductor C at zero potential, while the 
Coulomb system S lives on the line y = W {W > 0). 

We consider two cases : either the conductor C is an ideal one from the beginning, or 
it is the high-density limit of a Coulomb system (in that limit, the microscopic scale goes 
to zero). We compute the potential and charge correlations and the average stress tensor 
in each case. These quantities can be obtained by a macroscopic approach [6] , using linear 
response theory and a conducting behaviour assumption. Alternatively, exact microscopic 
results can be derived in a special model. 

A position r = {x, y) is conveniently represented by the complex number z = x -\- iy. 

2.2. The ideal conductor approach 

The line j/ = is assumed to be a grounded ideal conductor C, in that sense that the 
electric potential is constrained to vanish on that line. Thus, in the region |/, y' > 0, the 
potential at z' due to a unit point-charge at z is 



^(r, r') = - In 



z — z' 



z — z' 



(2.1) 



which does vanish on the hne y = 0. A Coulomb system S lies on the hne y = W. This 
model has been already studied [7]. Under the assumption that 5 is a conductor, a ma- 
croscopic approach [6] allows to compute the correlation function for the electric potential 
$ at two points r and r' (provided W is macroscopic and r and r' are at macroscopic 
distances from the line y = W). These correlations are (< ... means a truncated 
statistical average and /? is the inverse temperature) 



P < $(r)$(r') = - In 



z — z' sinh 2^ (-2 



z — z' sinh -^{z — z') 
if {y,y') e]0,W[^ , 



(2.2a) 



/3 < $(r)$(r') = - In 



z — z 



z — z' 



if (y,2/') G]0,W^[x]W^,+oo[ , 



(2.2b) 



/5 < $(r)$(r') = - In 



z — z 



iW 



z — z' 



if {y,y')e]W,+oo[' 



(2.2c) 



In the region j/ < 0, $(r) = without fluctuations. 

From (2.2), one can compute the correlations of the electric field E(r), and from the 
discontinuities of Ey{y) on the lines y = and y = W one obtains the correlations for the 
charge densities cr(r) (charge per unit length) on these lines : 



P < (J{r)a{r') > 



1 



27r2 



1 



2W J \smh^{x - x') 
if r and r' are on the ideal conductor, 



X — X' 



(2.3a) 



/3 < a(r)a(r') >t= - '- 



1 



2 \ 2W J \ cosh 2^7(3^ -a;') 
if r is on the ideal conductor and r' on the Coulomb system. 



(2.3b) 



and 



/3 < C7(r)c7(r') >■ 



1 



TT 



1 



2VFy \ sinh - x') 

if r and r' are on the Coulomb system. 



+ 



X — X' 



(2.3c) 



[(2.3b) and (2.3c) disregard the microscopic detail ; (2.3c) must be regularized x — x' = 
0]. As expected, if oo, the fluctuations (2.3a) on the ideal conductor C now alone in 

space disappear, as well as the correlations (2.3b) between the ideal conductor C and the 
Coulomb system S, while the fluctuations (2.3c) on the Coulomb system S become those 
of one conducting line [6] . 

From the averages < E^Ei, > associated to (2.2), one can compute the Maxwell stress 
tensor. Its only non-zero component is Tyy. At any point r between the lines y = and 
y = W, assuming that there is no potential difference between the lines, one found [7] 

TT 



yy 



(2.4) 



—Tyy is the force per unit length exerted on the Coulomb system S. 



2.3 The high-density limit approach 



We now assume that both lines y = and y = W are Coulomb systems, and the high- 
density limit is taken on the line y = 0. In that limit, the macroscopic approach [6] for 
computing the potential correlations is valid under the same conditions as in section 2.2 : r 
and r' should be at macroscopic distances from the line y = W, but there is no restriction 
about their distances to the line y = 0. One finds 

/3 < $(r)$(r') >^ = - ln\z- z'\ 

if {y,y') G] -oo,0[2 , 



(2.5a) 



P < $(r)$(r ) >- 



-In 



sinh 2|^(z - z') 



if iy,y') e]0,W[' 



(2.5b) 



/3 < $(r)$(r') = - In 1^ - z'\ 

if {y,y') e]0,W[x]W,+oo[ 
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(2.5c) 



P < $(r)$(r') >^ = - In 1^ - ^' - iW\ 



if {y,y') e]W,+oo[^ , 



(2.5d) 



P < $(r)$(r') >^ = -ln\z- z'\ 



if (2/,2/0 e] -oo,0[x]W,+oo[ . 



(2.5e) 



The difference between the electric potential correlations in this case and the previous 
case [equations (2.2)] is the electric potential correlation of a system where there is only- 
one Coulomb system on the line y = [6]. This difference raises from neglecting the 
fluctuations of the conductor at y = in the ideal conductor case. 
From (2.5) we obtain the charge correlations 



P < a{r)a{r') >^ = 



27r2 



TT 

2W 



sinh 2^(x — x') 



+ 



if y = y' = o , 



X — X 



(2.6a) 



P < a{r)a{r') >^ = -- ( — V ( ^ 

^ ^ ^ ^ 2 \2WJ l^cosh^O 



2iy{x a;') 
if y = and y' = W , 



(2.6b) 



P < a{r)a{r') >^ 



27r2 

if y ^ y' 



TT 

w 



sinh 2|7(a; - x') 



+ 



X — X' 



(2.6c) 



[eqs. (2.6) disregard the microscopic detail ; (2.6a) and (2.6b) must be regularized at 
x-x' = 0]. 
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The charge correlations in the Coulomb system at y = W are the same in both ap- 
proaches [equations (2.3c) and (2.6c)] as it was noticed before in other models [4, 2]. Also 
the correlation between a point on the conductor [y = 0) and a point on the Coulomb 
system (y = W) is the same in both approaches [equations (2.3b) and (2.6b)]. But the 
correlation between two points on the conductor at y = differ when the conductor is 
ideal [equation (2.3a)] and when it is the high-density limit of a Coulomb system [equa- 
tion (2.6a)]. 

The stress tensor is the same as in equation (2.4). This shows that we have the same 
results in both approaches for the force exerted on the Coulomb system S. What is special 
to the present model and what is general will be discussed in section 3. 

2.4- A solvable model 

The charge correlations (2.3) and (2.6) can be checked on a solvable microscopic model. 
In this section we consider that the Coulomb system 5 is a one-component plasma : the 
system is composed of particles of charge q moving in a rigid charged background. The 
two-dimensional one-component plasma is a solvable model in several geometries [8-11] 
when (3q^ = 2. 

a) The system such that C is an ideal conductor has been solved in [12, 13, 7], in 
the grand canonical ensemble. Let —qr] be the background charge density of S and ( the 
fugacity. The number density n and charge correlation in the Coulomb system S are given 
in terms of 




,ikx 



(2.7) 



as 



n = g{0) 



(2.8) 



and 



< a(r)a(r') >^= -q^\gix - x')f + q\S{x) 



(2.9) 
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It has been shown in [7] that, in the macroscopic hmit {r]W >> 1), these results agree 
with those from the macroscopic approach of section 2.2 ; (2.9) becomes (2.3a). 

b) The system such that C is a conducting Coulomb system can also be solved exactly. 
Now each line is a one-component plasma. The background charge densities are —qr/o for 
the line y = and —qij for the line y = W. Here we work in the canonical ensemble. Let 
N be the total number of particles. We consider first that we have two concentric circles, 
with radii R and R + W, on which each plasma lies, and then take the limit R ^ oo, with 
N = 27rr]oR + 2TTri{R + W), which ensures the overall neutrality. Adapting [11] by treating 
the radial coordinate r as a discrete variable which can have the values R and R + W, we 
introduce the N orthogonal functions 

M^) = {a 6r,R + Sr,r+w) , OKiKN-l , (2.10) 

where z = r e'^^ , 6 is the Kronecker symbol, and a is a (positive) parameter which controls 
[14] how the N particles are distributed between the two lines. The density n and charge 
correlations are given in terms of the projector 

Pir, r') . E (2.11) 

1=0 J2roe{R,R+w} lo \'^e{ro,9o)\^rod9o 

as 

n(r) = P(r,r) (2.12) 

and 

< a{r)a{r') >^= -q^\P{r,r')f + q\{r)6ryS{r{e - 9')) . (2.13) 

In the limit R ^ oo the two circles become two parallel lines. In this limit it is useful to 
define k = £/R. A summation over £ becomes an integral over k times R. We change our 
system of coordinates: let x — R9 and y — r — R. The projector becomes 

P(r, r) =Sy^y'{6y^oPi{x - x') + Sy^wP2{x - x')) 

+ {Sy,oSy',w + Sy,wSy',o)P3{x — x') (2.14) 



with 



27r(r?o+J?) 



dk 



P2{X 



1 + a-2 62^^= 



1 + a^e-^w^'^ 



]^ /•27r(77o+77) 



27tJo 



dk 



(2.15a) 



(2.15b) 



(2.15c) 



For a comparison with previous results, we define an alternative control parameter ^ by 



27rC = a-^e^'^^"^ 



(2.16) 



and keep Q fixed as we vary the other parameters 770, r/, W . 

When the density rjo of the conductor C becomes infinite, from (2.15b) where we make 
the change of variable k 277(7/0 + vi) ~ ^1 we obtain 



\P2{X)\ 



-,—ikx 



dk 



(2.17) 



27rio 1 + (27rC)-i e^w^C^-^'^^) 

Since |-P2(a;)| as given by (2.17) is identical to \g{x)\ as given by (2.7), the density and the 
charge correlation function on the Coulomb system S are indeed identical whenever C is 
an ideal conductor or the high-density limit of a Coulomb system. 

A more detailed comparison can be made when W is macroscopic. Then, neglecting in 
(2.15) terms ex.p{—r]oW) and exp{—r]W), after an averaging over oscillations of microscopic 
wavelength we obtain 



\P,{X)\' ~ \P2{X)\'' 



1 



47r2 



+ 



TT 



TT 



2W^sinh^ 



47r2 \ 2W^cosh^ 



(2.18a) 



(2.18b) 



The charge correlations obtained from (2.18) agree with the macroscopic ones obtained in 
(2.6). 
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3. General case 

In this section we consider the general case in d > 2 dimensions. The conductor C 
has any shape and the Coulomb system S occupies some region of space outside C. The 
Coulomb potential is 



To start with, C itself is considered as a Coulomb system with internal degrees of freedom. 

We shall use several quantities related to the electric potential correlations for the con- 
ductor C alone in space. Let ^c{^) be the electric potential at r created by the conductor 
C alone and let < ... >q be a truncated statistical average computed with the Boltzmann 
weight of the conductor C alone. The correlation < $c(r)$c(r') >o can be computed by 
linear response [6] ; it is related to the average electric potential change at r when a unit 
charge is put at r'. This potential change can be computed by macroscopic electrostatics. 
If the conductor is grounded, there are two cases : 

a) If r (or r', or both) is (are) inside the conductor. 




(3.1) 



< $e(r)$c(rO >l= G'o(r-rO 



(3.2a) 



b) If r and r' are outside the conductor. 



/? < $e(r)$c(r') >J= G*(r,r') 



(3.2b) 



where G* is defined by 



Ar [Go(r - r') - ^^(r, r')] = -^id5{v - v') 



(3.3) 



for r and r' outside the conductor, with ^2 = 27r, /is = 47r, ... fid = {d — 2)2n^/'^ /T{d/2) 
if (i > 2, and the condition Go{r — r') — G*(r, r') = if r (or r') is on the surface of the 
conductor ; Go{r — r') — G*(r, r') is the electric potential at r created by a unit charge 
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at r' in presence of a grounded ideal conductor. The expressions (3.2) which disregard the 
microscopic detail become exact in the limit when the correlation length goes to zero. 

Another remark useful for the following sections is that the fluctuations of the electric 
potential of a conductor C are Gaussian [6, 15, 16]. 

Let R be the set of particle coordinates of the conductor C, {rj} the set of the 
N particle coordinates of the Coulomb system S, dV the element of phase space of the 
Coulomb system and Hq{R) the Hamiltonian of the conductor. The total energy of the 
system (S plus C) is 

N 

iy({ri},R) = ^g,$c(r„R)+ J] Q,g,G'o(r, - r,) + ifo(R) , (3.4) 

i=0 ^<i<j<N 

where and qi are the position and charge of the i^^ particle of S. There might also be 
some short-range interaction between the particles of S, but we do not explicitly write it 
in (3.4) just for having a simpler notation. 

In the following sections we shall compute the partition function, the correlations, 
the stress tensor and the force exerted on the Coulomb system S, in the limit when the 
charge correlation length of the conductor C goes to zero (we shall call this limit the good 
conductor limit) and compare the results to the ones when the conductor is ideal (i.e. the 
potential on it is fixed, say to zero, without fluctuations). For the sake of simplicity we 
shall only treat in detail this case of a grounded conductor, but similar results hold for an 
insulated conductor (see section 3.5). 

3. 1 Partition function and statistical averages 

The partition function of the total system can be written as 

Z = j dT< e-/3E,^^*c(r,) >^ e-/3E,<,^^^.-«o(r^-'-,)^^ (3 5^ 

where < ... >o means the average over R with the Boltzmann weight exp(— /3i?o(R)) and 
Zq = J exp{—PHo{TC))dTi is the partition function of the conductor C alone. 
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Now, since the fluctuations of $c are Gaussian, 



< exp 



-/3 J]gi$c(ri) 



>o = exp 



= exp 



(3.6) 



where we have used (3.2b). Thus, the partition function becomes 



Z = Z*Zn 



(3.7) 



where Z* is the partition function of the Coulomb system 5' in presence of an ideal con- 
ductor, 

= I dV e-^^«« (3.8) 



where 



N 



H^s{vu..,VN) = --Y,iG*{ri,Vi)+ 9i9i[G'o(ri-r,)-G'*(ri,r,)] . (3.9) 

l<i<j<N 



i=l 



Hgff is indeed the standard Hamiltonian used in the ideal conductor approach. For instance, 
in the case of a plane ideal conductor, G* is the particle-image interaction ; it should be 
noted that the interaction —q'^G*{ri, r^) of a particle with its own image carries a factor 
1/2 in (3.9). 

The total free energy is F — F* + Fq where F* is the free energy of S in the presence 
of an ideal conductor and Fq the free energy of the conductor C alone. This was noticed 
previously in [4] for the model of a two-dimensional plasma near a metallic wall. 

Let A({ri}) be a microscopic quantity that does not depend on R. Its thermodynamic 
average is 

<A>= — fdr< e-^E.<^^^c(r.) ^-PE.,,1^'l.Go(r.-r,)^ 

Z* J 

— ^ /" (^Pg-/3iJeff(ri,..,riv)^ 
Z* 



-< A >efr 



(3.10) 
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where we have used (3.6). Thus the average of A can be computed by assuming from the 
beginning that the conductor C is ideal. 

3.2 Electric potential correlations 

Equation (3.10) does not apply to the electric potential correlations because the mi- 
croscopic electric potential is different in the cases of a good conductor or an ideal con- 
ductor. For the good conductor case, the microscopic electric potential is 




(3.11) 



while for the ideal conductor case it is 




(3.12) 



The average electric potential in the good conductor case is 




+ <$c(r)e-^S.'?^*c(rO 



>o e 



} . (3.13) 



Since the fluctuations of $c are Gaussian, 



< $c(r)e 



>o= -P^qi< $c(r)$c(ri) >o e 



^1^^ <li<lj<^c{ri)^cirj)>o 



(3.14) 



Using the covariance (3.2), one finds, for all r. 



< $(r) >=< $id(r) >eff 



(3.15) 



where < ^id(r) >eff is the average electric potential in the ideal conductor case. 
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We can compute the electric potential correlations in the same way : the correlation 
function in the good conductor case can be written as 



P < $(r)$(r') >= ^ / dT{ ^ q,q^Go{T - r,)Go{r' - rj)e- 

i 

+ J2qiGo{r' - r,) < $e(r)e-^S. >o 

i 

+ < $c(r)$c(r')e"^ >o ] } • (3.16) 

Using (3.15), 

< $e(r)$c(r')e-^^i^^*=^'''^ >o=[p^J2mj < >o< ^c(r')*c(r,) >o 

+ < $c(r)$c(r') >o ] 

[also a consequence of $c being Gaussian], and the covariance (3.2), we find 

< $(r)$(r') >'^= < $id(r)$id(r') >^^^+ < $c(r)$c(r') >J . (3.18) 

Thus, the correlation function in presence of a good conductor is the correlation func- 
tion in presence of an ideal conductor plus the correlation function for the good conductor 
alone in space. 

This is what was noticed in the example of section 2. 

3.2 Charge correlations 

If we are interested in charge correlations in the Coulomb system S, equation (3.10) 
applies because the microscopic charge density outside the conductor 

N 

p{r) = J2QA^-^^) (3-19) 
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does not depend on the coordinates R ; thus the charge correlations inside S are the same 
in both approaches. The surface charge density on C is given by the discontinuity of the 
normal electric field, thus using (3.15) we find that the average charge density is the same 
in both approaches. The same holds for the correlation between the density in S and 
the surface charge density on C. But, if we are interested in charge correlations on the 
conductor C, the correlations are different in the two approaches. Using equation (3.18) 
we can compute the difference in the electric field correlations and from it the difference in 
the charge correlations on the surface of C ; this difference is the surface charge correlation 
on C when it is alone in space. 

The example of section 2 illustrates these general results. 

3.4 The stress tensor and the forces exerted on the Coulomb system 

The Maxwell stress tensor is 

T^, = fi^' < E^E, - -^B^ > (3.20) 

where E = — is the electric field. Let V be some volume outside the conductor C. 
The total average electric force on V is 

F = / T.d'^-^S = / V.T(r)d^r . (3.21) 

JdV Jv 

It can be shown that this force is the same in both models although the electric 
potential correlations are different and consequently the stress tensor might be different. 
Indeed, from equation (3.18), the difference between the stress tensor in the good conductor 
case and the ideal conductor case is 

- 7;t(r) = l^d' - ^'Y^'^K) G*{r, r) (3.22) 

where d (resp. d') means partial differentiation with respect to the first (resp. the second) 
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argument of G* . The difference of the divergences is 



V(r) - 7;t(r) 



(3.23) 



Now, from (3.3) 9^9^G'*(r, r) = 0, and since G* is symetrical 9^G*(r, r) = 9j^G*(r, r). 
Thus, (3.23) is zero and the force is the same in both models. 

In the example of section 2, the stress tensor itself was the same in presence of either a 
good conductor or an ideal conductor. But this was an effect of the very peculiar symmetry 
of the model (invariance by translations along the x axis). 

A direct evaluation of the average force exerted on one particle of the system confirms 
that the two approaches give the same result. Let c/rjv-i be the phase space element of 
the system S when the i^^ particle is fixed and Z'^_.^ the partition function of the system S 
in presence of an ideal conductor when the i^^ particle is fixed. The average force exerted 
on the i^^ particle in presence of a good conductor is 



^qmGo{ri - Tk) 

.Kk 



+ <Vi 



-.(3.24) 



In the last term 



= -V,. 



kl 



Y^^l QkQtG*{rk,ri) 



Using (3.25) in (3.24) gives 



< Fi >=< Fj.d >efl . 



(3.25) 
(3.26) 



We obtain the same average force by both approaches. 
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3.5 The insulated conductor case 

In the former calculations we assumed that the conductor C was grounded. The 
same calculations can be carried out if the conductor C is insulated. Equations (3.7), 
(3.10), (3.15), (3.26) still hold for the insulated conductor case. The expression of ffeff 
now is a different one but it is still the Hamiltonian of the system S in presence of an 
ideal conductor : everywhere G*(r, r') must be replaced by G*(r, r') + Q{r) V{r') where 
Q{r) is the charge created by influence of a unit charge at r on the grounded conductor 
and V{r') is the electric potential at r' created by the conductor carrying a unit charge. 
The important fact to notice is that the relation (3.18) between the different potential 
correlation functions is still valid. 

4. Conclusion 

Two different methods for treating the problem of a Coulomb system near a conductor 
have been compared. 

The first method, where the conductor is considered from the beginning as ideal, ne- 
glects all fluctuations in the conductor. The second method treats the conductor as a 
conducting Coulomb system the charge correlation length of which goes to zero. Even in 
that limit, the fluctuations in the conductor do not vanish. This modifles the electric po- 
tential correlations by a term given by the electric potential correlation when the conductor 
is alone, but the average electric potential is not modifled. Because of this, quantities such 
as the electric fleld and charge density will have the same properties : their average is the 
same in both methods but their correlations differ by the correlation for the conductor 
alone. However, in the case of the charge correlations, this difference vanishes outside the 
conductor. 

The free energy in the good conductor case is just the sum of the free energy in the 
ideal conductor case plus the free energy when the conductor is alone in space. 
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The average stress tensor is modified only by a term the divergence of which is zero, 
so the average force exerted on any part of the Coulomb system is not modified. However, 
the fiuctuations of these forces will in general be different in the two methods. 
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